When finite-size corrections vanish: 
The S —1/2 XXZ model and the Razumov-Stroganov state 



o 
o 

(N 

c 
In 



•4— > 

CO 

+-* 



I 

C 

O 

o 



> 

m 
O 

m 

O 

On 
O 



Leonardo Banchi, 1 ' 2 Filippo Colomo, 2 and Paola Verrucchi 3, lj 2 
1 Dipartimento di Fisica, Universita di Firenze, Via G. Sansone 1, 1-50019 Sesto Fiorentino (FI), Italy; 



2 Istituto Nazionale di Fisica Nucleare, Sezione di Firenze, 
Via G. Sansone 1, 1-50019 Sesto Fiorentino (FI), Italy; 
Centro di Ricerca e Sviluppo SMC dell'Istituto Nazionale di Fisica della Materia 
Sezione di Firenze, via G. Sansone 1, 1-50019 Sesto Fiorentino, Italy. 



CNR, 



We study the one-dimensional S = 1/2 XXZ model on a finite lattice at zero temperature, 
varying the exchange anisotropy A and the number of sites N of the lattice. Special emphasis is 
given to the model with A = 1/2 and N odd, whose ground state, the so-called Razumov-Stroganov 
state, has a peculiar structure and no finite-size corrections to the energy per site. We find that 
such model corresponds to a special point on the A-axis which separates the region where adding 
spin-pairs increases the energy per site from that where the longer the chain the lower the energy. 
Entanglement properties do not hold surprises for A = 1/2 and N odd. Finite-size corrections to 
the energy per site non trivially vanish also in the ferromagnetic A — > — 1 + isotropic limit, which 
is consequently addressed; in this case, peculiar features of some entanglement properties, due to 
the finite length of the chain and related with the change in the symmetry of the Hamiltonian, 
are evidenced and discussed. In both the above models the absence of finite-size corrections to the 
energy per site is related to a peculiar structure of the ground state, which has permitted us to 
provide new exact analytic expressions for some correlation functions. 



I. INTRODUCTION 

Low dimensional magnetic systems have been acknowl- 
edged as intriguing physical systems for decades and still 
attract much interest, both from the theoretical and the 
experimental point of view. Many reasons justify such 
interest, and one more has been recently added, namely 
the possibility to use S — 1/2 spin models as tools for 
studying problems related to quantum information the- 
ory and quantum computation [3, 0j [H . Amongst one- 
dimensional systems, a preminent role is played by the 
Heisenberg Hamiltonian 



H 



(JxSi s i+1 



(1) 



where i runs over the sites of a chain, and Si are angular 
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momentum operators satisfying [Sf, Sj] 



(a, /3,7 = x, y, z). Models described by the Hamiltonian 
(fTj) constitute a class which is characterized, even at zero 
temperature, by the possible occurrence of peculiar phe- 
nomena, such as quantum phase transitions [H, satura- 
tion [1, [f| 0] , or factorization [9[ . Whether a specific 
model displays one such phenomenon depends on the de- 
tails of the exchange interaction and, in case, on the value 
of an external magnetic field. The values of the exchange 
parameters J x , J y and J z define each model, which in fact 
may get its name after such values. 

In this paper we will refer to the S =1/2 XXZ chain 
with a finite number of sites N and periodic boundary 
conditions, whose Hamiltonian can be written in the di- 
mensionless form 



H 



N 



i=l 



, (2) 



where erf are the Pauli matrices for the spin sitting at 
site i. The above expression is obtained from Eq. (fl]) 
by setting J x = J y = —4 and J z = 4A. The choice of 
the minus sign in front of the exchange interaction on 
the xy-plane implies no loss of generality in the thermo- 
dynamic limit, due to the possibility of changing such 
sign at will via a unitary transformation. However, when 
finite chains are considered, special care is due to this 
aspect, as explained in Section ITT1 

The behaviour of the model depends on the value of 
the anisotropy parameter: isotropic ferro- (A = —1) or 
antiferromagnetic (A = 1); Ising-like (|A| > 1); critical 
(|A| < 1). Within each interval, the actual value of A 
is not of particular interest, at least as far as the gen- 
eral phenomenology is concerned. However, in the case 
of the XXZ model Eq. @ there exists an exception to 
this statement: In fact, for A = 1/2, periodic bound- 
ary conditions, and a finite and odd number of sites, 
the ground state of the model, often referred to as the 
Razumov-Stroganov state, shows very peculiar features 
[ToL [Til . [lH which still stand as an unintelligible occur- 
rence. Such features do not fall within the framework of 
quantum phase transitions, as the model lies well inside 
the |A[ < 1 region of the XXZ Hamiltonian; moreover, 
the matter is relevant only as far as the number of sites 
of the chain is finite, so that the difference between odd 
and even number of sites stays meaningful. 

The Razumov-Stroganov state has been studied by sev- 
eral authors (see for example Refs. fl3l [fH flU [l6| ) . with 
different approaches and in many different frameworks, 
but somehow neglecting the fact that the corresponding 
model is not isolated in the phase diagram of the XXZ 
model. 

In this paper, we tackle the problem from a different 
point of view: Given the fact that a very peculiar ground 
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state occurs only for A = 1/2 and odd number of sites, 
we develop a comparative analysis of the behaviour of 
other XXZ models, i.e. models defined by Eq. |J2) with 
A 7^ 1/2, and N both even and odd, looking for clues 
about what in fact makes the model whose ground state 
is the Razumov-Stroganov state so special. From a phys- 
ical perspective, one of the most interesting feature of the 
Razumov-Stroganov state is the absence of finite-size cor- 
rections to its energy per site. This property is shared 
by the fully separable ground state of the finite-length 
ferromagnetic isotropic (A = —1) chain and also by the 
ground state of the A — » — 1 + limit, whose structure is 
however not as trivial, and deserves special attention. 

In the above framework, we have analyzed the effects 
of the finite length of the chain on the energy, the cor- 
relation functions, and some entanglement properties. 
The analysis presented is based on exact (i.e. available 
with infinite precision) results as far as the Razumov- 
Stroganov state and the ground state of the A — > — 1 + 
model are concerned, and on numerical data in all the 
other cases. Numerical diagonalization techniques un- 
derlie all our outcomes. 

The structure of the paper is as follows: In Section HT1 
we introduce the Razumov-Stroganov state and briefly 
recall its structure and properties. Some new exact re- 
sults for the two-point correlation functions of the corre- 
sponding chain (A = 1/2 and N odd) are also presented. 
In Section IIIII we consider the energy and the correlation 
functions at T — of several XXZ models, with |A| < 1 
and N both even and odd, focusing our attention upon 
finite-size corrections. In Section IIVI the same type of 
analysis is proposed for some entanglement properties at 
T = 0. Conclusions are drawn in Section fVl 

Some exact analytic expressions for the correlation 
functions of the A = 1/2 and the A — > — 1 + models 
are presented in the two Appendices. 



II. SPECIFICITY OF THE VALUE A = 1/2 

Let us consider the XXZ Hamiltonian on a finite 
chain with periodic boundary conditions and N odd, at 
the particular value A = 1/2: 



H = 



N 

E 



i 



(Of <7f +1 + O-M+l) + ^K+l 

&n+i = o\ ; 



i+i, ■ 2 - 
N odd. 



(3) 



It is worth noticing that, due to the specific conditions 
on the lattice, the negative sign of the transverse cou- 
pling cannot be reversed at will: Indeed, the rotation of 
7r around the z-axis of all the spins sitting at every other 
site can be safely performed only in the case of N even, 
or in the thermodynamic limit. 

The ground state of the above Hamiltonian is doubly 
degenerate for any N finite and odd, and the two de- 
generate ground states are eigenstates of YliLi with 



eigenvalue S^ ot = ±1/2. The two ground states are ob- 
tained from each other reversing the component along 
the quantization axis of each spin, and their structure is 
that of the so-called Razumov-Stroganov state. In the 
following we shall always refer to the Sf ot = +1/2 case. 

The Razumov-Stroganov state shows several intriguing 
properties [l(| EL E2] : i) its energy per site reads exactly 
E/N = — I, with no finite-size corrections; ii) the coef- 
ficients of the ground state on the standard basis (i.e. 
the basis where all operators o?, i = 1, ... N are diago- 
nal), are integer multiples of the smallest one; Hi) some 
of these integer numbers have a non trivial combinatorial 
interpretation. 

A complete understanding of these features is still 
lacking. It is worth recalling that the XXZ chain with 
A = 1/2 is known to have a similarly special ground 
state in two other cases: a) twisted boundary conditions 
and N even [l(| E3; V) open boundary conditions for 
whatever N, even or odd tialigi. 

The fact that, when suitably normalized, the coeffi- 
cients defining the ground state of Eq. §5§ on the stan- 
dard basis are all integers (see point ii) above) allows 
for their exact numerical evaluation [l9l . |2C| , for N not 
too large, limited only by computing capabilities. In- 
deed we have computed numerically these integer-valued 
coefficients with infinite precision, on a standard desktop 
computer, for system sizes up to N — 25. From the exact 
knowledge of the ground state the correlation functions 
are readily determined and, as in the present case they 
necessarily assume rational values, we have computed 
them with infinite precision, again for chain lengths up to 
N = 25. We report in Appendix IA1 for illustrative pur- 
poses, some of our results for the second neighbour lon- 
gitudinal two-point correlation function (ofof +2 ) N . Here 
and below (...)„ denotes the expectation value over the 
ground state of a chain of N spins. Further results are 
available upon request [l||. 

It appears that the knowledge of the exact values of 
such correlation functions for a finite set of values of N, 
together with the fact that these are rational numbers, al- 
lows to determine their analytic expressions as functions 
of N (obviously, only valid for odd N). This is a rather 
exceptional situation for an interacting critical system. 
In particular, from the infinite-precision data mentioned 
above, and the exact thermodynamic-limit values com- 
puted in Refs. [2l| and [13], we have been able to work 
out exact analytic expressions for the r th -neighbour lon- 
gitudinal and transverse two-point correlation functions, 
for r = 1, 2, 3, 4, 5, and arbitrary N. The expressions for 
r = 1 coincide with those computed by means of purely 
analytic methods in Ref. [ll[ 



16N 2 ' 



(°f°f+l); 



2N 2 



(4) 

while those for r — 2,3,4,5 were previously unknown. 
We report them in Appendix lAl 

Let us briefly illustrate the method used in deriv- 
ing such expressions, taking as an example the case 
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of the second-neighbour longitudinal two-point correla- 
tion function (0faf, 2 ) N , whose exact values for N — 
3, 5, ... , 15, are reported in Appendix [A] Subtracting 
from such values the corresponding exact value in the 
thermodynamic limit (in this case, 7/64, see [2l|), and 
factoring into prime numbers the resulting denominators, 
it is easy to infer for these denominators a behaviour, 
as a function of N, of the form 2 6 A 2 (A 2 - 4). The 
corresponding numerators are also expected to behave 
as a polynomial in A, although of lower order, since 
the fraction as a whole should vanish in the thermo- 
dynamic limit. The first four values of (of <J? +2 ) N , for 
N = 3, 5, 7, 9, are sufficient to determine the polynomial 
in the numerator, thus completely fixing the expression 
of (crfal +2 ) N for arbitrary A, see Eq. (|A2]) . The ob- 
tained expression reproduces exactly the available nu- 
merical values of {<jI<jI +2 ) n1 f° r N larger than 9. From 
a strictly mathematical point of view, the results of such 
procedure can only be conjectural, though strongly sup- 
ported by the numerical data. 

The above procedure can be applied to r th -neighbour 
correlation functions, but the degree of the involved poly- 
nomials increases rapidly. For generic values of r, the de- 
nominators appears to behave as jV 2 [( r+1 )/ 2 l n!=i'(-^ 2 ~ 
4i 2 ) r ~ 2l+1 . When r = 6, for example, the denomina- 
tor is of degree 24 in A, and even under the reason- 
able assumption (supported by all lower-r examples) that 
the numerator is an even polynomial in A, we need 12 
data to determine it, and possibly one more to check 
the result. These data can be taken from the values for 
N = 7, 9, . . . 31, which however lie beyond our computing 
capabilities. In this respect it is worthwhile to empha- 
size that, although the derivation of the results presented 
above and in Appendix [X] rely heavily on computer aided 
evaluations, they are all exact. 

III. OTHER MODELS 

In this Section we aim at getting some deeper insight 
into the physical mechanisms possibly related to the pe- 
culiar ground state of Hamiltonian ([3]); in particular, we 
wish to highlight the specific role played by each of the 
two conditions defining the model, i.e. A = 1/2 and A 
finite and odd (in the case of periodic boundary condi- 
tions) ; to this purpose we develop a comparative analysis 
of the behaviour of other XXZ models, i.e. models de- 
fined by Eq. © with A ^ 1/2, for different values of 
A. 

We first recall that, despite the most spectacular fea- 
tures being observed for A odd, the model ([3]) has a 
precise specificity also in the thermodynamic limit. In 
fact, as demonstrated by Baxter in Ref. [23], whenever 
the exchange parameters of an infinite Heisenberg chain, 
Eq. (fT]), satisfy the condition 

JxJy "T" JxJz ~\~ JyJz = (5) 

the energy per site E/N of the ground state gets the 



value 

\{J X + Jy + J z ), (6) 

in the case of the XXZ model © , this may only occur for 
A = 1/2, a condition which therefore defines a somehow 
special point, though by no means related with quantum 
critical transitions, nor with possible factorization or sat- 
uration of the ground state. 

In Section [TT] we have seen that one of the most in- 
triguing peculiarities of the model Eq. © is that finite- 
size corrections do not enter the energy expression, which 
necessarily implies (H) N /N — —3/2, given that condition 
(|S"|) is fulfilled. Due to translation invariance, and given 
that TL = J2i T~U,i+i, the absence of finite-size corrections 
is a property of the local energy, by this meaning that the 
expectation value of the nearest-neighbour interaction is 
not affected by the possible addition of no matter how 
many spin-pairs (though adding one single spin would 
drastically change the whole picture, bridging the model 
to the essentially different case of A even); the above 
feature is seen to follow from the exact cancellation of 
the A-dependent terms in the nearest-neighbour correla- 
tion functions, when combined as necessary: In fact, by 
Eqs. flU), it is 

(Hi,i+i) N = ^ 4 ( cr i fcr i+i)jv + -^\ a i G l+\)N 
5 3 1 3 3 
~4 ~ 4A 2 ~~ 4 + 4A 2 ~ ~2 ' ^' 

which shows that, despite the correlation functions being 
affected by finite-size corrections, the interaction energy 
is not. 

Let us now analyze how XXZ models behave for dif- 
ferent values of A and A. For the sake of clarity, the 
finite-size corrections relative to any physical observable 
O will be hereafter studied in terms of the quantity 8 N (O) 
defined by 

5 N (0)=(0) N - lira {0) M . (8) 

We first consider S N (Wj^+i) which is shown in Fig. [1] for 
different values of A and as a function of A: The qual- 
itative difference between even (upper panel) and odd 
(lower panel) number of sites is evident: finite-size cor- 
rections are negative for all values of A and A, except 
A > 1/2 and A odd. 

In other terms, while in general the energy per site 
increases with the size of the chain, a special situation 
occurs for odd A and 1/2 < A < 1, where adding spin 
pairs lowers the energy: Within the odd- A sector, Hamil- 
tonian §3§ is hence found to correspond to the separating 
case between these two qualitatively different behaviours 
of the finite-size corrections in the energy of the ground 
state. The Insets in Fig. Q] show how E/N varies with 
N for various anisotropics in the vicinity of the value 
A = 1/2: It is evident that such value is crossed with 
continuity, both in the odd- and in the even-A sector; 
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FIG. 1: 5 N (Hi, i+ i) versus N for A = 0.7, 0.5, 0.3, 0, -0.3, 
-0.5, -0.7. The full line is a guide for the eyes and connects 
data for A = 1/2. Insets: E/N versus N for A= 0.48, 0.49, 
0.495, 0.499, 0.50, 0.501, 0.505, 0.51, 0.52. The anisotropy 
parameter A grows as indicated by the arrows. Upper panel: 
N even; lower panel: N odd. 



this particularly means, in the former case, that finite- 
size corrections smoothly vanish for A — > 1/2. A simple 
perturbative analysis for very small e = A — 1/2 yields 



E 
N 



3e 
2PP 



(9) 



in complete agreement with the above scenario. 

It is to be noticed that the existence of a unique value 
of A, where S N (7ij,i+i) changes its sign for all values of 
N odd, is definitely not granted, as shown below in the 
case of the nearest-neighbour correlation functions along 
the ^-direction. 

Before considering the correlation functions, we remind 
that our results are obtained via the numerical diagonal- 
ization of the Hamiltonian, so that not only some physical 
observables, but all the coefficients of the ground state 
on the standard basis are available, for N both even and 




N odd 




FIG. 2: Nearest-neighbour correlation functions {cf<J~ +1 ) N , 
and | (erf af + i) N | versus A for different values of N odd and 
even. N grows as indicated by the arrows, from 5 to 25. The 
Inset shows the difference (af cr% + i) N — (of of +1 )g as a function 
of A, for N=5, 7, 9, and 25. 



odd. This allows us to notice the following general fea- 
tures of the ground state of the XXZ model with | A| < 1: 
i) all the elements of the standard basis belonging to the 
sector with S* ot = 1/2 for N odd, or S* ot — for N even, 
enter the decomposition of the ground state; ii) the co- 
efficients on the standard basis grow with the number of 
antiparallel adjacent spin pairs featuring the correspond- 
ing element; Hi) keeping the minimum coefficient fixed 
to unity, all the other coefficients grow with A; iv) the 
coefficients are all equal in the A — » — 1 + limit. 

We further underline, due to its relevance for the 
following discussion, that the nearest-neighbour corre- 
lation functions along the z-direction are negative not 
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only for A > 0, as expected, but also for —1 < A < 
0. Indeed, although the exchange interaction in the z- 
direction is ferromagnetic, configurations where adjacent 
spins have opposite components along the z-axis are en- 
ergetically favoured by the predominant exchange inter- 
action in the rry-plane, which is (afaf +1 + o\a\x^) = 
2(a^~a^ +1 +a^<T^ +1 ). For the sake of clarity, in the lowest 
panels of Fig. Owe plot the absolute values of (cf(rf +1 ) Jv . 

We now focus upon the finite-size corrections to the 
quantities (^fcr^ + i) N and (crfa? +1 ) N : From Fig. [2] it is 
clear that {ufa^ +1 ) N for odd N behaves in a peculiar way: 
in fact, while S N {^t a 7+i) i s always positive (the longer 
the chain the less correlated adjacent spins are on the 
xy-plane), we find that for any odd N it exists a value of 
A where S N (afaf +1 ) changes its sign, so that, as far as 
the anisotropy is larger than such value, adjacent spins 
get more and more correlated along the z-direction as the 
length of the chain grows, while the opposite occurs oth- 
erwise. It is to be noticed that, at variance with the case 
of 5 N (Hi t i+i), the value of A where S N (a? of +1 ) vanishes 
does depend on N. This numerical observation is made 
evident by plotting the difference (of of +1 ) N — (of af +1 ) N , 
for different values of N and fixed N': should finite-size 
corrections vanish at a precise value of A indipendently 
on N, all the lines in the Inset of Fig. [2] (where we have 
chosen N' = 9) would cross at the same point, which is 
seen not to be the case. 

Let us now discuss the above results. We have found 
that nearest-neighbour correlation functions along the z- 
direction in the S — 1/2 XXZ model with A > — 1 arc 
always negative due to prevailing role of the exchange 
interaction on the xy-plane, which evidently favours con- 
figurations with antiparallel adjacent spins. On the other 
hand, the necessary occurrence, in the case of N odd, 
of at least one pair of adjacent spins with the same 
component along the z-direction, locally frustrates this 
propensity to antiparallelism. As such constraint con- 
cerns just one spin pair, no matter the length of the 
chain, it may only affect finite-size corrections, and its 
relevance is smeared out as N increases. This mecha- 
nism of local frustration is responsible for the anomalous 
behaviour of finite-size corrections to (ofof+i}^ for N 
odd, and specifically for the change of sign of S N (of of +1 ), 
which ultimately makes it possible for 5 N {TLi t i+i) = 
-4Ma+ar+i) + AM«+i) to vanish at A = 1/2, 
and stay positive for A > 1/2. The model whose ground 
state has the specific structure predicted by Razumov 
and Stroganov is thus found to be a special point in the 
finite-TV XXZ phase-diagram, as it separates the region 
where adding spin pairs increases the energy per site from 
that where the longer the chain the lower the energy per 
site. 

From this point of view, the XXZ model in the A — > 
— 1 + limit is similarly special (for a precise definition of 
such model, see Appendix iBj) . Indeed, its ground state 
has a peculiar structure: It is a superposition with iden- 
tical coefficients of all the elements of the standard basis 
within the S£ ot = 1/2 (Sf ot = 0) sector, for N odd (even). 



This allows, as shown in Appendix[Bl for the exact evalu- 
ation of the correlation functions. In particular, similarly 
to Eq. ([7]), one find, for example for N odd, 



-1 



1 
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1 

N 



-1 



(10) 



and analogously for N even. We thus see that finite- 
size corrections to the energy per site (but not to other 
physical quantities) do vanish for this model as well. 

We notice that there exists another class of models 
whose expression for the energy per site is given by 
Eq. ([5]), namely any antiferromagnetic XYZ Heisenberg 
system on a bipartite lattice with periodic boundary con- 
ditions, in the presence of a uniform magnetic field whose 
value fulfills the condition given by Kurman et al. in 
Ref. Q. When such condition is fullfilled, the ground 
state factorizes, meaning that it gets the unexpected 
classical-like structure of a tensor product of single-spin 
states. This phenomenon has recently attracted much 
attention, due to its relevance as far as the entanglement 
properties are concerned 0, [U, [2f| [2(| . In this respect, 
one should notice that while the absence of finite-size cor- 
rections in the models studied by Kurmann et al., as well 
as in the A = — 1 ferromagnetic isotropic model, directly 
follows from the factorized structure of the ground state, 
thus characterizing the expectation value of whatever 
physical observable, in the Razumov-Stroganov state, as 
well as in the A — > — 1 + ground state, finite-size correc- 
tions vanish only in the energy per site, following a pre- 
cise cancellation in the expectation value of the Hamil- 
tonian. 



IV. ENTANGLEMENT PROPERTIES 

The investigation of entanglement properties has re- 
vealed a powerful tool in studying spin models, par- 
ticularly in unveiling peculiar properties of the ground 
state ll HI |H Sfl HI, HI , as in the case of the XYZ 
antifcrromagnet in the uniform field mentioned above [9|] . 
Several authors have specifically addressed the analysis 
of the XXZ model in terms of entanglement properties 
[H, HH, [H, H3, [H, H^, and more recently the case of 
A = 1/2 and N odd has been studied in such context 
[2fj| . In this framework, the idea that the peculiar struc- 
ture of the Razumov-Stroganov state might be related 
with properties of some entanglement measure is sug- 
gestive. We have therefore developed the same type of 
analysis presented in Section IIII1 for the bipartite entan- 
glement between two 5=1/2 spins separated by r lattice 
spacings, as measured by the concurrence [4(| El, H3| 



C r>N = 2max{0, C' N } , 



(11) 



where 



N 



(12) 
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FIG. 3: ft versus A and N 



We have numerically evaluated C r! jv for various A, r, 
and N. In Figs. [3] and 0] we show C\,n and C2,Ar versus N 
and A. Apart from very small oscillations between odd 
and even N, essentially due to different finite-size cor- 
rections in the correlation functions and to the fact that 
(of ) 2 N in Eq. (JT2J) is positive for N odd and vanishes for 
N even, no peculiar behaviour is observed. The concur- 
rence between nearest neighbours C\^m is always finite 
and is larger for larger A, i.e. for more marked antifer- 
romagnetic exchange along the z-direction (as expected 
after the analysis presented in Ref. [ZD, HH). All pairs 
of non-adjacent spins are disentangled in the Razumov- 
Stroganov state, no matter the value of N, but this is not 
a specific feature of the model (J3]), as already shown by 
the overall picture presented in Ref. 3l|. After all, and 
despite the particular structure of its ground state, the 
case A = 1/2 and N odd is not found to display peculiar 
features, at least as far as pairwise entanglement is con- 
cerned. The way transverse and longitudinal correlation 
functions combine in Eq. (|12p does not cause any relevant 
change in the finite-size corrections of the concurrence. 

Let us now consider the other point where finite-size 
corrections to the energy per site vanish, i.e. the ferro- 
magnetic isotropic point A = — 1. As a matter of fact, 
the model in the A — > — 1 + limit displays features which 
deserve some comments. In the thermodynamic limit, all 
the C r are expected (Hf to switch on as A — > —1+ which 
implies, due to the monogamy inequality [IfJ, |4Q , that 
all the concurrences vanish at the ferromagnetic isotropic 
point A = — 1, consistently with the fact that the ground 
state of the XXZ chain with A < — 1 is fully separable. 
On the other hand, when N is finite (no matter if even 
or odd), not only all the C rtN switch on as A — > — 1+, 
but they all flow to the same value, as seen in Fig. In 
particular, from the exact evaluation of the correlation 




FIG. 4: C 2 versus A and N 




-1.0 -0.5 0.0 0.5 -1.0 -0.5 0.0 0.5 



A A 

FIG. 5: C' riN as a function of A for N = 25 (left) and N = 9 
(right) for r < 5. 



functions reported in Appendix [B] one can find that 




(13) 

Moreover, when A < — 1 the result C t .n = for all r 
stays valid also for finite N; therefore, and at variance 
with the thermodynamic limit, a discontinuity of all the 
concurrences occurs at A = — 1. 

We underline that the one-tangle [45[ r = l-Eol 11 ")™ 
in the XXZ model with periodic boundary conditions and 
|A| < 1 is different from zero and does not depend on A 
(being r = 1 — l/N 2 for N odd, and r = 1 for N even), 
while for A < — 1, r vanishes. The observed discontinu- 
ity at A = — 1 (at variance with the one characterizing 
the concurrence) is not related with the finite size of the 
system, as it obviously survives in the thermodynamic 
limit. It is rather connected with the onset of critical 
behaviour for |A| < 1, which implies the vanishing of 
the magnetizations and a consequently maximal entan- 
glement content. 
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We understand the mechanism leading to such dis- 
continuities as follows: From our results we observe, as 
stated in Section IIII| that the ground state of the XXZ 
model for A — > — 1 + may be written as P + Hi\ip)i where 
P + is the projector over the Hilbert subspace of the chain 
corresponding to Sf ot — 1/2, and \cp) is a single-spin pure 
state which may be chosen at will, provided it is not one 
of the two eigenvectors of o~ z (so as to ensure the pro- 
jection does not vanish). The notation indicates that all 
the spins along the chain are in the same single-spin pure 
state: In fact, the fully separable state n,|y>)j is just one 
of the infinite ground states of the ferromagnetic XXX 
chain. We therefore see that it is the projection over the 
5 t z ot = 1/2 Hilbert subspace which injects entanglement 
into the fully separable ground state of the ferromagnetic 
XXX model. We notice that such projection is not a lo- 
cal operation, and it can hence modify the entanglement 
content of the system. The above reasoning is straight- 
forwardly extended to the TV even case, replacing P + 
with P°, i.e. with the projector over the S% ot = Hilbert 
subspace. In other terms, the easy-plane character of the 
XXZ model with |A| < 1, which is embodied in the con- 
dition S% ot = 1/2 for TV odd, or 5f ot = for TV even, 
causes the ground state to develop long-ranged pairwise 
entanglement as the ferromagnetic isotropic point is ap- 
proached from above. 

Since pairwise entanglement between two spins does 
not exhaust the possible types of entanglement, and there 
might be other measures that exhibit a special behaviour 
near or at A = 1/2, we have numerically evaluated the 
entanglement entropy for a wide range of bipartition of 
the full chain, and for various TV and A. In the par- 
ticular case where the bipartition consists of two blocks 
of adjacent spins, we find a very good agreement with 
the theoretical predictions of conformal field theory (4?J . 
Such agreement is observed not only for A = 1/2 (as al- 
ready emphasized in Ref. [2(|) but all over the interval 
— 1 < A < 1, and TV both even and odd. It is worth 
underlining that such agreement is observed to hold for 
blocks of at least two spins with great accuracy, even for 
rather short chains, with TV as small as 5. 



V. CONCLUSIONS 

Our study of the ground state structure as related with 
finite-size effects in the one-dimensional 5 = 1/2 XXZ 
model with periodic boundary conditions in its critical 
phase (T = and |A| < 1 ) has highlighted some rele- 
vant features; these particularly emerge for A = 1/2 and 
TV odd, as well as in the ferromagnetic isotropic limit, 
A — > — 1 + (for TV both even and odd). Indeed, these two 
models share a common property, namely the fact that 
finite-size corrections to the energy per site vanish, de- 
spite being finite as far as other physical quantities are 
concerned. 

In the former case, we have found that the model rep- 
resents a special point in the TV-odd sector (where its 



ground state has the structure of the Razumov-Stroganov 
state), as it separates the region where adding spin pairs 
increases the energy per site of the chain, from that where 
the reverse holds. We think this result, besides giv- 
ing some further insight into the behaviour of the mag- 
netic XXZ model, could get a specific meaning for the 
fermionic models whose Hamiltonian is mapped into the 
5 = 1/2 XXZ one (see for example Ref. 14]), possibily 
in terms of a physical mechanism related with a properly 
defined chemical potential. 

In the A — > — 1 + limit, the peculiarities of the model 
do also arise from the very specific (despite not as intrigu- 
ing as that of the Razumov-Stroganov state) structure of 
the ground state, which guarantees the exact cancella- 
tion of finite-size corrections to the energy per site, for 
TV both even and odd. However, in this case the peculiar 
ground state structure most relevantly reflects on bipar- 
tite entanglement: all the concurrences are finite and get 
the same value, no matter the distance between the two 
spins considered. In fact, we have found that there is 
a region, while approaching the ferromagnetic isotropic 
point from above, where any two spins along the chain are 
entangled, given that the finite size of the chain prevents 
the monogamy of the entanglement to force the vanishing 
of all the concurrences; in such region the entanglement 
between a selected spin pair can be varied not only by 
tuning A but also by properly choosing TV, a possibility 
which might be relevant in the design of quantum devices. 

Finally, we recall that the two above mentioned mod- 
els, besides sharing a somewhat specific finite-size be- 
haviour, present another rather exceptional peculiarity, 
namely the possibility of accessing the exact structure of 
their ground state. Indeed, in both cases, the (suitably 
normalized) coefficients of the ground state in the stan- 
dard basis are all integers: This naturally sets the models 
into the realm of combinatorics, and has permitted us to 
derive new exact analytic expressions for their correlation 
functions. 
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APPENDIX A: SOME EXACT RESULTS AT 

A = 1/2 

We report here the exact values of the second 
neighbour longitudinal two-point correlation function 
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(« +2 >^ for TV = 3, 5,..., 15: 



(°f°f+2)3 = 


1 

~3 


( fT i cr l+2)5 = 


1 

25 


K^2>7 = 


4 

49 


( cr i :cr l+2)9 = 


28 
297 


( CT i CT ?+2>ll = 


157 
1573 


( CT i CT ?+2>13 = 


191 
1859 


(^^+2)15 = 


1732 
16575 



(Al) 



Note that all these values are rational numbers. This re- 
markable property allows for their infinite-precision nu- 
merical evaluation, and permits to derive the exact ex- 
pressions of the correlation functions for arbitrary (odd) 
values of A. We report below our results for the r th - 
ncighbour longitudinal and transverse two-point correla- 
tion functions, r = 2, 3, 4, 5: 



«^2>, 



( cr i fcr t+2 )iv = 



7 



41 

28 



3 / 227 + 227V 2 



26 \N 2 (N 2 -4) 



105 
"28" 



1 + 2N 2 
A 2 (TV 2 - 4) 



(A2) 



(A3) 



401 45 /-2205 + 53247V 2 + 267V 4 + 2127V 6 \ 
K^i^i+zlN - ~2}2 + 212 y TV 4 (TV 2 - 4) 2 / 



(A4) 



+ 4399 9 

\ a i a i+3/N - -^15- + ^15 



-11025 + 12767V 2 - 73587V 4 + 45167V 6 
TV 4 (TV 2 - 4) 2 



(A5) 



184453 



2 22 2 22 TV 4 (TV 2 - 4) 3 (TV 2 - 16) V 



( 16094658825 - 4071808726TV 2 



+ 1675416103TV 4 - 269157300TV 6 + 41497687TV 8 + 2281766TV 10 ) (A6) 



, , _ . 1751531 



224 224^4 (^2 _ 4)3 ^ 2 _ ^ V 



4695690825 + 2413434266TV 2 



95214949 



2 31 



+ 



231TV 6 (TV 2 - 4) 4 (A 2 - 16) 2 V 



276248393A 4 + 556663980A 6 - 134362697A 8 + 10926614A 10 ) (A7) 



( 3019990307709375 - 1994398549102575A 2 



+ 22018044125468TV 4 - 149133792747084A 6 + 28928020131522A 8 - 4476210910162A 10 

+ 693645483372TV 12 - 43191843804TV 14 + 2051652263TV 16 ) (A8) 



3213760345 



■ (4227986430793125 - 565146993503925TV 2 



2 35 2 35 A 6 (A 2 - 4) 4 (TV 2 - 16) 2 

- 383849970492812TV 4 - 183487079646084A 6 + 70690548132342A 8 - 35899621486502TV 10 

+ 8073828981732A 12 - 796229443764A 14 + 28577677213TV 16 ) (A9) 
I 



APPENDIX B: THE A -1+ MODEL 

We provide here some analytic results on the XXZ 
chain of length A, in the A — > — 1 + limit. It is worth em- 



phasizing that, by construction, this model differs from 
the isotropic A = — 1 ferromagnct. It is defined as a limit 
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from the critical region, and its ground state, together 
with the expectation value of all physical quantities, is 
also to be intended as the result of this limit. In partic- 



tot 



1/2 for N odd, 
-1 + e, there ob- 



ular, the model is critical, with S 1 
and Sf 0t = for TV even. Setting A = 
viously are finite corrections to all the exact expressions 
we provide below, but they all vanish as e — > 0. 

It is easily seen, from analytic considerations, and nu- 
merical experiments, that the ground state of the XXZ 
chain in the A — > — 1+ limit is a superposition with iden- 
tical coefficients (which we set to unity) of all the ele- 
ments of the standard basis corresponding to S^ ot = 1/2 
for TV = 2n + 1 (odd), and to S£ ot = for TV = 2n (even). 
Such elements are herafter indicated by \ei). Their num- 
ber is (^Q ; the norm of the state consequently reads 

TV = \J (n)- Due to the peculiar structure of the ground 
state, the evaluation of correlation functions reduces to 
a simple combinatorial enumeration, as shown below. 

Let us consider the correlation functions on the xy- 
plane: it is (<J+aJ) N = J\f~ 2 ^2i m (ei\a^ aj \e m ) , and the 
elements \e m ) that do contribute to the double sum are 



just those with af\e m ) 



-\e m ), and a z Ae m ) = \e m ); 



each of them will contribute exactly 1 and their number 
is (n-i) i so that 



TV + 1 

4TV ' 
N 

4(TV - 1) 



TV odd, (Bl) 
TV even, (B2) 



As for the correlation functions along the z-axis, it is 
( a i a j) = N~ 2 J2i m ( e i\ a i (7 j\ e m)i and each term in the 
sum is easily seen to contributes with ±<5; m . The various 
terms can be collected, according to the orientation of 
the i th and j th spins in each \e m ), into four groups, so as 
to obtain 

w)„ = at 2 r- 2 ) - ( n -_\) ( N n z\) + r„- 2 )] 

(TV - 2n) 2 - TV 



TV(TV-l) ' 



(B3) 



and hence 



TV odd, 



(B4) 



l i'n TV — 1 ' 



TV even, (B5) 



no matter the distance between site i and j. 



no matter the distance between site i and j. 

The fact that two-point correlation functions do not 
depend on the distance is a direct consequence of the 
particular structure of the ground state, whose non van- 
ishing coefficients on the standard basis are all identical. 
It is worth mentioning that each of these correlation func- 
tions assumes the same value on the chain of length 2n — 1 
and 2n. Finally, notice that these correlation functions 
have different values with respect to the A = — 1 isotropic 
ferromagnetic model. The consequent discontinuity van- 
ishes in the thermodynamic limit. 
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